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The achievement of half-metaUicity with ferromagnetic (FM) coupUng has become a key technology for the 
development of one- dimensional (ID) nanoribbons for spintronic applications. Unfortunately, in previous 
studies, such a half-metallicity always occurs upon certain external constraints. Here we, for the first time, 
demonstrate, via density functional theory (DFT), that the recent experimentally realized gallium sulfide 
nanoribbons (GaSNRs) can display an intrinsic half-metallic character with FM coupling, raised from 
Ga-4s, Ga-4p and S-3p states at the Ga-dominated edge. Furthermore, the novel half-metallic behavior with 
FM coupling here is rather robust, especially for GaSNRs with large width and thickness, and can be 
sustained to the room temperature. Thus, our results accidentally disclose a new ID spin nanomaterial, 
which allows us to go beyond the current scope limited to the graphene, boron nitride (BN), zinc oxide 
(ZnO) and molybdenum sulfide (M0S2) nanoribbons, toward more realistic spintronic applications. 

Since the discovery of ID nanoribbons, it has sparked an intense research effect toward the understanding of 
such materials with the promising application in nano spintronic devices \ Nevertheless, such proposed 
application requires the ability of 100% spin-polarized at the Fermi leveP. A half-metallicity can filter the 
current into a single spin channel, and thus fully meets this demand^'^. Consequently, the achievement of half- 
metallicity becomes the forefront of spintronics study. 

To date, various strategies have been proposed theoretically to realize such ID half-metallic nanostructure. Son 
et al. predicted the half-metallicity in zigzag graphene nanoribbon (GNR) under a sufficient transverse electric 
field across the ribbon^ This may open an exciting pathway to next-generation spintronics. However, in contrast 
to the widely applied gate voltage or bias electric fields in the current technology, the required transverse electric 
field for the half-metallicity transition increases with increasing ribbon width, and thus its realization is challen- 
ging^. To circumvent this difficulty, hydrogenation has also been proposed to realize the half-metallicity the- 
oretically in GNR^, while the experimental manipulation of hydrogen atoms to precise positions is still a challenge 
because hydrogen atoms are easy to form in a random way on a host surface. Recently, edge modification receives 
a lot of attentions due to the effective realization of half-metallicity in ID nanostructures^"^l Although theoretical 
method seems successful in showing half-metallicity, it remains experimentally impractical because of the 
difficult in current edge micromanufacturing techniques. Another alternative approach to realize half-metallicity 
is to apply external strain. For example, the half-metallicity in M0S2 nanoribbon (M0S2NR) can occur when a 
tensile strain in zigzag direction arrives 11%^^. However, the experimental achievement still represents a difficulty 
because the electronic and magnetic properties of M0S2NR exhibit a very sensitive response to the applied strain 
that needs the fine control of strain on M0S2NR. In fact, other various routes, such as injecting charges^^'^^, doping 
atoms or organic molecules^^'^^, inducing nanoholes^^'^^, and hybriding different materials with the similar 
structure^'^^"^^ also have been proposed to realize the novel half-metallicity in nanostructure, but an intrinsic 
half-metallicity without any external constraints is not yet demonstrated and is more desirable for the practical 
spintronic applications. In this respect, the effective acquisition of a ID material with intrinsic half-metallic and 
FM character represents a key challenge to realize its practical application. 

As a kind of ID nanostructure, GaSNRs have been experimentally fabricated^^'^^. Shen et al. successfully synthe- 
sized high-quality ID GaSNRs via a simple vapor- solid method, and pointed out that the morphology and structure 
of the products can be easily controlled by substrate temperature and evaporation source^^ By a simple thermal 
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Figure 1 | Top and side views of the optimized geometric structure of 
Z^^iQj^i-GaSNR. The ribbon is periodic along the x direction, and its 
width is denoted by w in the y direction. The region contained by the 
dashed rectangle is a unit cell used to calculate the magnetic moment, 
energy difference between AFM and FM couplings, and electronic 
properties. 

evaporation process, the well-crystalline ultrathin GaSNRs have also 
been experimentally realized on silicon substrates, and are proposed as 
a promising electron source^^. As a consequence, the experimental 
realization of GaSNRs might resolve the key challenge mentioned 
above. However, the scientific issue is if such GaSNRs can exhibit 
the intrinsic half-metallic character that has not been found in pre- 
vious ID nanostructures. If such a half- metallic property indeed exists, 
the effect from ribbon's width and thickness, especially at the room 
temperature, must be fully explored for practical applications. In order 
to unlock their valuable properties, a thorough understanding of struc- 
tural, electronic, and spin properties of GaSNRs is highly desirable. 
However, to our best knowledge, the ID GaSNRs with unique prop- 
erty have not yet received any theoretical attention so far. 

In this paper, we explored the electronic and magnetic properties 
in a series of zigzag GaSNRs using DFT calculations, and found that 
all GaSNRs exhibit intrinsic half-metallic character. Moreover, the 
induced half-metallicity presents a substantial FM character through 
long-rang magnetic coupling, especially for GaSNRs with large width 
and thickness. Hence, our results provide a new pathway to explore 



GaSNRs-based spintronics. In contrast to previous studies, our 
results exhibit the following advantages over the existing ones: (1) 
it is not necessary to apply an external electric field or strain to 
GaSNRs, through which the half-metallicity can be found; (2) it is 
not necessary to induce point or line defects in GaSNRs, around 
which the half-metallicity may appear; (3) it is not necessary to 
modify the zigzag edge of GaSNRs by doping foreign atoms or 
organic molecules, near which the half-metallicity may be induced. 
Although these methods can be achieved experimentally, until now 
the effective manipulability at nanoscale is still a challenge. 

Results 

Geometries of zigzag GaSNRs. Our systematic study began with the 
geometric detail of a single gallium sulfide (GaS) layer. The perfect 
GaS layer is fourply with two Ga atoms in the middle which are 
covalently bonded to S atoms located in top and bottom sheets, 
forming S-Ga-Ga-S arrangement^^'^^. Computed lattice constant of 
3.58 A, S-Ga distance of 2.37 A, and Ga-Ga distance of 2.46 A are in 
good agreement with previous results^^'^^. The zigzag GaS 
nanoribbon (GaSNR) here was constructed by cutting the GaS 
layer along two parallel zigzag lines. For convenience, a zigzag 
GaSNR with a width of "n" and layers of "m" is described as 
Z^=n,i=m-GaSNR. Figure 1 shows the top and side views of a fully 
relaxed Z^=io, i=i-GaSNR (as an example). It is apparent that the 
honeycomb network and fourply structure in Z^=io,i=i -GaSNR are 
well remained, supporting the experimental realization of such a 
GaSNR^^'^^. Previous studies have suggested that such unsaturated 
ID nanoribbons are closely associated with magnetism^ Indeed, 
our calculations demonstrate that Z^=io,i=i -GaSNR can lead to the 
magnetic moment of 0.753 |iB. Here the unpaired spins of Z^= io,i= i" 
GaSNR are mainly contributed by Ga (0.240 iib) and S (0.120 iib) 
atoms at the Ga-dominated edge, as shown in Figure IS in the 
supplementary information. In contrast, no magnetism was found 
in all armchair GaSNRs, and thus we only consider zigzag GaSNRs in 
the following study of magnetic behaviors. 

Half-metallicity of zigzag GaSNRs. One of the most interesting 
applications associated with the magnetism in nanoribbon is the 
novel half-metallicity. In order to explore if such a zigzag GaSNR 
can possess the novel half-metallic character in detail, we plot 
Z^=io,i=i-GaSNR's spin-polarized band structure and density of 
states in Figure 2(a) and 2(b), respectively. Indeed a unique half- 
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Figure 2 | (a) Spin-polarizedband structure of Z^=io,i=i-GaSNR, (b) spin-polarized total density of states (TDOS) of Z^=io,i=i-GaSNR and projected 
density of states (PDOS) of Ga and S atoms and (c) spin-unpolarized band structure of Z^=io,i=i -GaSNR. The Fermi level is set to zero. 
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Figure 3 | Top and side views of spatial spin distribution of 10,1=1- 
GaSNR (isovalue 2 X 10"' e/A'). 

metallicity is observed. The spin-down electrons are metallic with 
two levels (a and P) crossing each other at the Fermi level while the 
spin-up electrons are semiconducting due to the existence of a band 
gap as large as 1.0 eV. The high spin polarization at the Fermi level 
can ensure a 100% passage of preferred spin, and it opens a transport 
channel only for spin-down electrons and blocks the channel for 
spin-up electrons. Thus, the charge transport is mainly dominated 
by the spin- down electrons. Considering the effect of edge 
passivation on the half-metallicity, we tested the passivation of Ga 
atoms at the Ga-dominated edge by H atoms. With H passivation, a 
and P levels in spin-down band can be kept well and the half- 
metallicity is nearly unchanged. The interesting result can be 
attributed to the composition of a and p levels. Here both a and p 
levels are composed not only by Ga-4s and Ga-4p but also by S-3p 
states. Thus, when Ga atoms were passivated by H atoms, a and P 
levels can still be kept well due to the existence of S-3p states. It is 
different to the result found in zigzag BN nanoribbons (BNNRs) 
where the half- metallic character can be somewhat affected by the 
passivation of N atoms at the N-dominated edge by H atoms^\ In 
contrast, although previous studies also have predicated the half- 
metallicity upon certain external constraints mentioned above, 
such a half-metallicity without any external constraints here is 
intrinsic and thus should be more suitable for spin filter device 
applications. Considering the origin of special half-metallic 



character, we also calculated the spin-unpolarized band structure 
of Z.. = io,i=i"GaSNR in Figure 2(c). It can be seen that two half- 
filled states (HFSs) near the Fermi level were found. These HFSs 
near the Fermi level are zero-energy states, and the existence of 
degenerate HFSs is not stable. When the degree of freedom for 
spin is unrestricted, these half-filled HFSs split, and Z^= 10,1=1 - 
GaSNR exhibits a half-metallic character. As a matter of fact, the 
achievement of half-metallicity due to half-filled HFSs has been 
demonstrated in previous study^^. 

FM coupling of zigzag GaSNRs. To visualize the spin distribution of 
Z^= 10,1=1 -GaSNR, we impose its spin density, i.e. the charge density 
difference between spin-up and spin-down channels, as shown in 
Figure 3. We can see that the unpaired spin almost concentrates on 
the Ga and S atoms in the Ga-terminated edge while the contribution 
from the inner and S-terminated edge atoms can be almost negligible 
agreeing with the magnetic moment's analysis in Figure IS in the 
supplementary information. This is different to the finding in zigzag 
GNR and M0S2NR where the unpaired spins are contributed from 
both edges^^'^^, but similar to the result of zigzag ZnO nanoribbon 
(ZnONR) in which the unpaired spin is only contributed from the O- 
terminated edge^^. The origin of unique spin distribution can be 
attributed to the Ga-S interaction. After relaxation, we found that 
the Ga-S bond length in the S-terminated edge was reduced by 
0.15 A, which results in the enhancement of the Ga-S covalent 
bonding interaction and thus quenches the magnetism of Ga and S 
atoms. Conversely, in the Ga-terminated edge, Ga-S bond length was 
slightly increased by 0.05 A, resulting in the reduction of the Ga-S 
covalent bonding interaction and leading to the unpaired electrons 
accumulated on Ga and S atoms. In fact, the change of the magnetism 
due to the variation of covalent interaction is analogous to the result 
in previous studies^^'^^. Note that the unpaired spin electrons here 
align parallel, i.e. FM coupling, and exhibit a substantial collective 
character by long-range magnetic coupling. Considering the stability 
of such FM coupling, the energy difference between antiferro- 
magnetic (AFM) and FM couplings needs to be further 
investigated. Figure 2S(a) in the supplementary information gives 
ten possible AFM couplings, i.e. AFMj (1 < j < 10), for monolayer 
GaSNRs. Calculated energy differences between AFM and FM 
couplings are about 18, 19, 17, 16, 18, 19, 20, 18, 18 and 20 meV 
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Figure 4 | (a) Front view of spatial spin distribution of Z^=n,i=i -GaSNR (5 < n < 12) (isovalue 2 X 10"^ e/A^) (b) spin-polarized band structure of 
Zw= 12,1=1 -GaSNR. Note that the spin-polarized band structure of Zw=n,i=i-GaSNR (5 < n < 11) presents the similar half-metallic character that was not 
given. The Fermi level is set to zero. 



SCIENTIFIC REPORTS | 4 : 5773 | DOI: 1 0.1 038/srep05773 



3 



Table I | Magnetic moments of FM and AFM couplings and energy difference between ARM and FM couplings per unit cell for a series of 
Zw=n,i=rGaSNRs (5 < n < 1 2, 20). Here Mfm and Mafmi (1 — j — 1 0) denote magnetic moment of FM and AFMj couplings, respectively. 
AEj (1 < j < 1 0) denotes the energy difference between AFMj (1 < j < 1 0) and FM couplings 
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for AFMl, AFM2, AFM3, AFM4, AFM5, AFM6, AFM7, AFM8, 
AFM9 and AFM 10, respectively. Here calculated energy dif- 
ferences between AFM and FM couplings are based on a unit cell 
as shown in Figure 1. These results suggest that FM coupling is more 
stable than AFM coupling for all cases considered. Conventionally, 
the induced magnetism is often located on adsorbed atoms or near 
vacant sites. However, for a doped system, one of the crucial issues 
for its applications in spintronic devices is whether the local spin 
moments induced by defect states can lead to a collective magnetism, 
which is an essential requirement for any spintronic applications. 
Unfortunately, this critical issue, as well as its implications, was often 
overlooked in the previous studies of the doping-induced magnetism 
of nanostructures. Here the finding of FM spin coupling with a 
substantial collective character renders GaSNR a more practical 
candidate for spintronic applications. 



Influence of ribbon width on half-metallicity and FM coupling. 

Since experimentally synthesized nanoribbons always have a relative 
large width, we investigated the width effect on the electronic and 
magnetic properties on zigzag GaSNR. Interestingly, from our 
calculations, zigzag GaSNR with a large width still exhibits the 
half- metallic character with FM coupling. The relevant spatial spin 
distributions and the spin -polarized band structure are given in 
Figure 4(a) and 4(b), respectively. The interesting result can be 
attributed to the stable edge structure which dominates the half- 
metallic character and is independent of the GaSNR width. In 
order to further understand coupling behaviors, we calculated 
magnetic moments of FM and AFM couplings and energy 
difference between AFM and FM couplings for GaSNR with 
different width, as shown in Table I. It can be seen that FM 
coupling is more stable than AFM coupling for all cases 
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Figure 5 | (a) Front view of spatial spin distribution of Z^=5j=in-GaSNR (1 < m < 8) (isovalue 2 X 10"^ e/A^) (b) spin-polarized band structure of 
Z^=5,i=8-GaSNR. Note that the spin-polarized band structure of Zw=5,i=m-GaSNR (1 < m < 7) presents the similar half- metallic character that was not 
given. The Fermi level is set to zero. 
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Table II | Magnetic moments of FM and AFM couplings and energy difference between ARM and FM couplings per unit cell for a series of 
Zw=5 l=rn"GaSNRs(l < m < 8). HereMpMcmdMAFMi (1 — j — 1 0) denote magnetic moment of FM and AFMj couplings, respectively. AEj (1 
< j < 1 0) denotes the energy difference between AFMj (1 < j < 1 0) and FM couplings 
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considered. For FM coupling, Z^=5,i=i-GaSNR has a lowest 
magnetic moment with the relatively low energy difference, and 
Z^_2o,i=i-GaSNR has a largest magnetic moment with the 
relatively large energy difference. It is expected that magnetic 
moment and energy difference increase with the ribbon width, 
indicating that not only the magnetic moment of FM coupling but 
also the stability of FM coupling are enhanced with the increase of the 
ribbon width. However, no magnetism is found from the calculation 
of infinitely single-layered GaS because the magnetic moment per 
GaS molecular formula decreases with the increase of width, which is 
analogous to that studied previously in MoSaNR^^ Nevertheless, the 
experimentally synthesized GaSNRs with large width are expected to 
possess very large magnetic moment and stability for FM coupling 
and thus are expected for practical applications. 

Influence of ribbon thickness on half-metallicity and FM coupl- 
ing. Since the experimentally synthesized nanoribbons generally 
consist of several layers, we also considered the thickness effect on 
the electronic and magnetic properties of zigzag GaSNRs. Here the 
width of GaSNR was fixed to be w = 5 for computational expediency 
and a multilayer ribbon was constructed by placing the planar layers 
on top of each other with an AB stacking. As noted in Figure 5(a) and 
5(b), the spin density of Z^=5j=ni-GaSNR can be viewed as a simple 
combination of those in Z^= 5 1=1 -GaSNRs, and even with more 
layers zigzag GaSNR still presents a half-metallic character. As a 
comparison, zigzag ZnONR exhibits a different behavior, where 
the edge atoms of the top layer can covalently bind with the edge 
atoms of the bottom layer and thus the magnetism diminishes when 
the number of layers becomes even^^. This interesting result may be 
attributed to the very weak van der Waals interaction between 
GaSNRs where the average separation between two neighboring 
GaSNRs is about 3.4 A. In order to further understand coupling 
behaviors, we calculated magnetic moments of FM and AFM 
couplings and energy difference between AFM and FM couplings 
for GaSNR with different thickness, as shown in Table II. Here ten 
possible AFM couplings, i.e. AFMj (1 < j < 10), for multilayer 
GaSNRs are shown in Figure 2S(b) in the supplementary 
information. It can be seen that FM coupling are more stable than 
AFM coupling for all cases considered. For FM coupling, calculated 
magnetic moment for Z^=5j=in-GaSNR structure can be 



approximately determined to be 0.703 m jiB, and calculated energy 
difference between AFM and FM couplings for Z^=5^i=ni-GaSNR 
structure can be approximately evaluated to be 18 m, 19 m, 17 m, 
16 m, 18 m, 19 m, 20 m, 18 m, 18 m and 20 m meV for AFMl, 
AFM2, AFM3, AFM4, AFM5, AFM6, AFM7, AFM8, AFM9 and 
AFM 10, respectively. It indicates that the experimentally 
synthesized GaSNRs with large number of layers are expected to 
possess very large magnetic moment and stability for FM coupling 
and thus are expected for practical applications. 

Influence of temperature on half-metallicity and FM coupling. 

Considering whether the experimentally fabricated GaSNRs can be 
served at room temperature, spin-polarized molecular dynamics 
(MD) simulations were carried out for 1 fs at 300 K. Figure 6 
depicts the fluctuation of the total energy (Etot) and the magnetic 
moment (M) of three representative GaSNRs, i.e., Z^=5i=i-, 
Z^^ 12,1-1- and Z^=5 i=8-GaSNRs, as a function of time. We found 
that, after 2000 steps, the geometries of three configurations still 
remain, but with small variations in total energy and magnetic 
moment. This can be rationalized by the fact that the Ga-S binding 
energies in GaSNRs are much larger than the thermal energy 
corresponding to the room temperature. As a consequence, 
GaSNRs still exhibit half-metallic character with FM coupling. 
(Note that half-metallic character and FM coupling here were 
determined by DFT simulations at 0 K with the configuration 
obtained from MD simulations at room temperature.) Hence, the 
novel half-metallic character, along with FM coupling, renders the 
GaSNRs an excellent candidate for spintronic appHcations at room 
temperature. 

Influence of calculated functional on half-metallicity and FM 
coupling. For the exchange-correlation functional, DFT calculated 
band gaps are generally lower than experimentally determined band 
gaps. Hence, we have tested band structures in a series of zigzag 
GaSNRs with Heyd-Scuseria-Ernzerhof (HSE06) hybrid functional. 
Although HSE06 calculated band structures are somewhat different 
from those by DFT calculation, HSE06 calculation also gives the 
similar magnetic behaviors, i.e. half-metallicity and FM coupling. 
Hence, our calculated results survive to the choice of functional. 
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Figure 6 | The variation of the total energies, Etot> and the magnetic moment, M, with time from molecular dynamics simulations for Z^=5j= 
^=12,1=1- and Z^=5j=8-GaSNRs. 



Discussion 

Half-metallicity with FM coupling has become a key technology for 
ID nanostructure's application in spintronics. Previously, such a 
half-metallicity always occurs under certain external constraints, 
such as applying an external electric field or strain, inducing point 
or line defects, as well as modifying the zigzag edge. Nonetheless, 
considering the practical spintronic applications, the effective 
acquisition of an intrinsic half-metallicity without any external con- 
straints is particularly important, which represents a significant chal- 
lenge for modern technologies. Based on spin-polarized DFT, we, for 
the first time, demonstrated the existence of intrinsic half-metallicity 
with FM coupling in recent experimentally realized GaSNRs. 
Although corresponding AFM couplings present metallic character 
as shown in Figure 3S in the supplementary information, the unique 
half- metallic character with FM coupling here is more stable than 
metallic character with AFM couplings, especially for experimentally 
realized GaSNRs with large width and thickness. Furthermore, the 
novel half-metallic behavior with FM coupling can be sustained to 
the room temperature. Thus, our results accidentally disclose a new 
ID spin nanomaterial toward the more realistic spintronic applica- 
tions. In contrast to previous studies, our results exhibit the following 
advantages over the existing ones: (1) we do not need to apply an 
external electric field or strain in GaSNRs to realize their half-metal- 
licity; (2) we do not need to induce point or line defects in GaSNRs to 
realize their half-metallicity; (3) also, we do not need to modify the 
zigzag edge of GaSNRs by foreign atoms or organic molecules to 
realize their half-metallicity. Although these methods can been 
achieved experimentally, until now the effective manipulability at 
nanoscale is still a challenge. 

Methods 

The first principle periodic calculations based on DFT are performed using the 
Vienna ab initio simulation package (VASP). Within generalized gradient approxi- 
mation (GGA), we considered Perdew— Burke — Ernzerhof (PBE) exchange and 
correlation functional. The pseudopotentials with 4s^4p^ and 3s^3p^ valence electron 
configurations are used for Ga and S atoms, respectively. The electronic wave func- 
tions were expanded using a plane-wave basis set with a cutoff energy of 350 eV. The 
Brillouin zone was sampled by 20 special k points for atomic structure relaxation and 
30 k points for the electronic structure calculation. AID periodic boundary condition 
was applied along the ribbon. Both of vacuum distances between ribbon edges and 
between layers of two adjacent images are set to be about 10 A. Full optimization of 
the atomic structures including the atomic positions and lattice parameters has been 
carried until all the force components are smaller than 0.01 eV/A, and the conver- 
gence of the electronic self- consistent energy is less than 10"^ eV. DFT calculated 



band gaps are generally lower than experimentally determined band gaps, and thus 
HSE06 hybrid functional also has been used to verify the correctness of electronic and 
magnetic properties of zigzag GaSNRs. HSE06's results show that major magnetic 
behaviors, i.e. half-metallicity and FM coupling in zigzag GaSNRs, are similar to those 
found from DFT calculation. 
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